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We introduce a class of PT-symmetric systems which include mutually matched nonlinear loss 
and gain (in other words, a class of PT-invariant Hamiltonians in which both the harmonic and 
anharmonic parts are non-Hermitian). For a basic system in the form of a dimer, symmetric and 
asymmetric eigenstates, including multistable ones, are found analytically. We demonstrate that, if 
coupled to a linear chain, such a nonlinear 'PT-symmetric dimer generates new types of nonlinear 
resonances, with the completely suppressed or greatly amplified transmission, as well as a regime 
similar to the electromagnetically-induced transparency (EIT). The implementation of the systems 
is possible in various media admitting controllable linear and nonlinear amplification of waves. 
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Introduction. In the past few years, the study of 
systems exhibiting the parity-time (PT) symmetry has 
drawn a great deal of attention. The underlying idea 
is to extend canonical quantum mechanics by introduc- 
ing a class of non-Hermitian Hamiltonians which may, 
nevertheless, exhibit entirely real eigenvalue spectra pQ. 
A necessary condition for the Hamiltonian to be "PT- 
symmetric is that its linear-potential part V(x), being 
complex, is subject to the spatial-symmetry constraint, 
V(x) — V*(—x). The complex PT-symmetric potentials 
can be realized in the most straightforward way in op- 
tics, by combining the spatial modulation of the refrac- 
tive index with properly placed gain and loss [2]. This 
possibility has stimulated extensive theoretical [3] and 
experimental : 4\ studies. 

In the PT-symmetric Hamiltonians introduced in the 
context of the field theory and optics, solely the harmonic 
part features the matched gain and loss, while the anhar- 
monic part, if any, is Hermitian, giving rise to nonlinear 
dynamical models in which only the linear part features 
the balanced dissipation and amplification [5]. In this 
work, we put forward an extension of the PT symmetry, 
in the form of Hamiltonians whose anharmonic part too 
includes the mutually matched loss and gain. Solving the 
corresponding dynamical equations, we demonstrate that 
such nonlinearly-VT-symmetric systems also give rise to 
eigenstates with real frequencies. New findings, which 
are specific to the systems with the matched nonlinear 
gain and loss, are eigenstates with a spontaneously bro- 
ken spatial symmetry (that, as mentioned above, imple- 
ments the PT symmetry of the Hamiltonian) , and mul- 
tistability of eigenstates. A straightforward application 
of these states is realized by coupling the PT system to 
linear chains: we demonstrate that this setting gives rise 
to new types of multistable nonlinear Fano resonances, as 
well as to transmission regimes with a very strong ampli- 
fication, and those similar to EIT (electromagnetically- 
induced transparency) . Systems of this type can be im- 



plemented in optics, using saturable absorbers [5] and 
two-photon losses to realize the nonlinear PT symmetry, 
as well as in any medium which allows nonlinear ampli- 
fication of waves, such as cavity polaritons [7], surface 
plasmons jS], and magnons [5]. 

We start by introducing a solvable nonlinearly-PT- 
symmetric system, in the form of a dimer, in which the 
symmetric linear gain and loss terms come along with 
their nonlinear mutually conjugate counterparts, 

iip A = (E + ij - ij2\^A\ 2 + x\i>A\ 2 )i>A + Vip B , 

(1) 

iip B = {E- i'fo ~ ij2\ipA\ 2 + x\^b\ 2 )^b + VipA- 

Here the overdot stands for the time derivative, 70 > 
accounts for the linear gain and loss, which act on com- 
plex variables tpA and ipB , respectively, E is a frequency 
shift (with respect to a linear chain if the dimer is cou- 
pled to it, see below), 72 accounts for the PT-symmetric 
nonlinear loss and gain (as shown below, stable eigen- 
states are obtained with 72 > 0, i.e., if the nonlinear 
loss competes with the linear gain, and vice versa), \ is 
the strength of the nonlinear frequency shift, and V is a 
coupling coefficient. 

Symmetric modes. Symmetric eigenstates, with \4>a\ = 
\iPb\, are sought for as ipA,B(t) = Aexp (— iuot ± id/2), 
with the amplitude and phase shift determined by the 
following equations: 

[ x A 2 -(uj-E)] 2 + ( 1q - 12 A 2 ) 2 = V 2 , (2) 

tan 5 = ( 7o - 72 A 2 ) (uj-E-xA 2 )' 1 . (3) 

Depending on the parameters, Eq. ([2]) may yield no phys- 
ical solutions with A 2 > 0, a single solution (monostabil- 
ity), and bistability, with two physical roots. The bista- 
bility occurs under conditions 

(u>-E) 2 > V 2 - 1 l 

7072X" 1 + (u-E)> ^{uj-E) 2 +1 l-V\ 
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while the monostability condition is (ui — E) < V 2 — 
7q. Although the system is dissipative, the symmetric 
eigenstates with the real frequencies form a continuous 
family parameterized by arbitrary frequency u, which is 
a manifestation of the VT symmetry. 

Asymmetric states and multistability. The system ad- 
mits solutions with broken symmetry too (A ^= B): 
ij) AtB {t) = {Ae i5 / 2 ,Be~ iS / 2 } e~ luJ \ with 8 determined by 
the same equation ^ as above. Unlike the symmetric 
eigenstates, the asymmetric ones exist at a single fre- 
quency, which is typical to dissipative systems: 
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FIG. 1. (Color online) The linear chain with the side-coupled 
elements featuring the nonlinear VT symmetry. The arrows 
indicate incident, reflected and transmitted waves. 



was = E + (70/72) X- 
Further, amplitude A is determined by equation 



(5) 



(A 2 ) 2 - ( lQ / l2 )A 2 + V 2 (x 2 + J 2 )- 1 = 0, (6) 

cf. Eq. ([2]), while the other amplitude is given by 
B 2 = (70/72) — A 2 , which shows that the asymmetric 
eigenmode exists only for 72 > 0, i.e., when the nonlin- 
ear 'PT-symmetric loss/gain terms compete with their 
linear counterparts. We stress that the asymmetric solu- 
tions are supported by the balance of the nonlinear gain 
and loss, as they do not exist for 72 = 0. 

The above relations yield two physical solutions (i.e., 
the bistability) for the asymmetric modes, with A 2 , B 2 > 
0, at x 2 Il2 > 4 (V 2 /7q) — 1, and no solutions in the op- 
posite case. Under this condition, inequalities Q hold 
too for uj = was 1 he., the system gives rise to the multi- 
stability, with four coexisting eigenstates, two symmetric 
and two asymmetric. 

The scattering problem. The next step is to couple the 
dimer to a chain transmitting linear discrete waves ip n (£), 
as shown in Fig. [I] (note that the entire system remains 
PT-symmetric). Here we focus on the most fundamental 
version of the system, with x = 0, the nonlinearity be- 
ing represented by the matched cubic loss and gain, the 
respective coupled system being 

iip A = E^ A + i( 1() -j 2 \^ A \ 2 )^ A + Vipo, (7) 

ilpn = C (lpn-1 + 1pn+l) + V5 n>0 (lp A + 1p B ) , (8) 

iip B = Eip B - i (70 - "f2\ipB | 2 ) 4>B + VlpQ, (9) 

where C is the coupling constant in the linear chain. The 
general solution corresponding to the scattering of inci- 
dent waves with amplitude / on the VT complex is looked 
for as 



Vv, 



Ie i(kn-ut) + Re -i{kn+wt) ( n < 0) . 
Te i(kn-ujt) ( n > 0) , 



(10) 



where wavenumber k > is determined by the dispersion 
equation for the linear chain, k = cos -1 (lu/2C), while 
R and T are the amplitudes of reflected and transmitted 
waves. A straightforward analysis of Eqs. ( 8|10 ) at n — 



yields R = ipo — I , T — ipo, and the expression for tpa in 
do) . 

J A,B- 



terms / and tp^ 



ip =I + iV {2Csinky 1 (^ 0) + ^ 0) ) . (11) 



The substitution of expression (11) into the stationary 
version of Eqs. ^ and ^ leads to a system of complex 
cubic equations: 



(E - u) ?/4 0) b + W 2 (2C sin ky 1 ^ + ip 



±M 7o - 72 



Va.b 



A = -vi. 



One should solve Eq. (1 1 2h for ifi A B at given / and u> 



(12) 



Then, tp can be found trom Eq. (11), and, eventually, 



the reflection and transmission coefficients can be found. 

The scattering in the symmetric regime. In the lin- 
ear system (72 = 0), Eq. (12) yields only symmetric 
solutions, with \ip A \ — \iPb\- The corresponding scat- 
tering spectrum, displayed in Fig. [2j demonstrates two 
noteworthy effects. One is the suppression of the trans- 
mission by the degenerate side-coupled elements without 
the gain and loss, 70 = 0. In this case, the eigenfrequen- 
cies of both elements are identical, and their excitation 
results in the resonant reflection at uj = E, which can 
be explained in terms of the Fano resonance jlOJ . The 
presence of the weak linear gain and loss, with 70 <C 1, 
lifts the degeneracy between the attached sites, leading 
to a response resembling the electromagnetically-induced 
transparency (EIT) [11], with the total transmissivity 
(T = 1) at uj = E, between resonant reflections on the 
pair of slightly detuned linear VT elements. 

In the system combining the linear chain and the non- 
linear VT scatterer with 72 > 0, one can find symmetric 



solutions with ip 



(0) 



where d> is real. 



First, we consider the case of (j) ^ \J 70 ~pfi (this value 
plays a special value, as shown below). Then, the sym- 
metric mode exists at u) = E, with 4> determined by 
equation 



72 3 - 7o0 + VI = 0, 



(13) 
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FIG. 2. (Color online) The transmission coefficient for several 
values of the gain/loss factor 70 in the linear system (72 = 0). 
Other parameters are E — 0.1, V = 0.2, and (7 = 1. 



which yields a single real root for P m = I 2 > 
(4/27)7o/ (^ 2 72)j and three real solutions (tristability) 
in the opposite case. According to Eq. (11), all these 



solutions realize the perfect EIT-like transmissivity, with 
T = 1 [the horizontal blue line in Fig. |3ja)]. The family 
of the symmetric states is displayed by the blue curves 
in Fig. [3](b) , where the tristability occurs at Pi n < 1/27. 

Nonlinear Fano resonances. In contrast to its linear 
counterpart, the nonlinear system may support complete 
suppression of the transmission (T — 0), i.e., nonlinear 



Fano resonances [ID] . From Eq. \ll\\t follows that ip^ - 



V^ 0) = 2UCV- 1 sinfc for ip Q = T = 0. 



The substitution 



of this into Eq. ( 12 ) leads to the system 



(o) 

A.B 



±i 



7o 



;(0) 
72 IpA B 



Va.e 



0, (14) 



which gives rise to a continuous family of symmetric non- 
linear Fano resonances, with uj — E and 



^ { ab = ^70/72 exp (±iS) , 



cos S = (2V)- 1 v/(4C 2 - E 2 ) (72/70) P h 



(15) 



(16) 



(recall the above EIT-like symmetric family, correspond- 
ing to T = 1, had (f> 7^ \/7o772 and cos<5 = 0). Equation 
(16 1 imposes condition cos 2 5 < 1, i.e., 



Pin<W 2 (70/72) (4C 2 -E 2 y 



(17) 



Thus, at to — E, the continuous family of the symmetric 
nonlinear Fano resonances exists in this interval of the 
intensity of the incident wave. The novelty of the result is 
that the Fano resonance is usually obtained as an isolated 
solution. 



The asymmetric scattering regimes. Equation ( 14 1 for 
the nonlinear Fano resonances also admits two ultimate 
asymmetric states, with the vanishing excitation at one 
of the VT elements: lu = E and 



2,^ 



0, 



(18) 




FIG. 3. (Color online) The transmission coefficient (a) and 
excitation intensity of the side-coupled VT elements (b) for 
70 = 0.01, 72 = o'.OOOl, V = 0.2, C = 1, and u = E = 0.1. 
The nonlinear Fano resonances correspond to T — 0. The 
red and blue curves depict, respectively, the full multitude of 
asymmetric scattering regimes, produced by Eq. ( 12 1, and the 
symmetric one corresponding to Eq. (13 1. The inset in (b) 



zooms the tristability region in the latter case. The horizontal 
dotted line corresponds to Eq. ( 15 1. 



or vice versa, 
exists at R n = 



with A <^ B. In either case, this solution 



V 2 (70/72) (4C 2 



E 2 



This point falls 



into the range (17) of the existence of the symmetric 
Fano-resonance solutions, which implies intrinsic bista- 
bility of the nonlinear Fano resonances. 



Equation ( 12 1 gives rise to other asymmetric scattering 
regimes, which, in particular, may produce a strong res- 
onant amplification of the transmitted wave. The com- 
plete set of the asymmetric scattering states is depicted 
by the red curves in Fig. [3] 

Stability. The stability of the above analytical solu- 
tions was checked in direct simulations of Eqs. ([7|-(|9|. 
The results demonstrate that the ultimate asymmetric 



state ( 18 ) with the nonzero excitation at the linear-loss 
element is stable, while its counterpart with the excita- 
tion at the linear-gain element is unstable, transforming 
itself into an oscillatory mode (apparently, a limit cycle) , 
see left panels in Fig. [4] These results can be readily 
understood following the similarity to previously studied 
systems composed of coupled cores with the linear gain 
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FIG. 4. (Color online) The perturbed evolution of those asym- 
metric (left) and symmetric (right) Fano-resonance modes 
from Fig. [3] which are dynamically unstable. Filled circles 
indicate the initial states. 



and loss acting separately in them, which also give rise 
to a pair of stable and unstable modes [12] . 

Symmetric Fano-resonance modes ( 15 ) are unstable, 



also developing intrinsic oscillations, with a very low 
transmissivity and spontaneously broken symmetry be- 
tween the VT elements, \i))a\ ^ IV'fll) as shown in the 
right panels in Fig.[4j In fact, such asymmetric states cor- 
respond to the nearly perfect Fano resonance in Fig. [3] at 
Pi„ = 4. The resonantly-amplified transmission regimes 
are also unstable, due to enhancement of the field on the 
linear-gain element. 

Nonpropagating modes. The dispersion relation for 
Eq. ([8j demonstrates that frequencies of the propagating 
waves belong to the respective phonon band, |w| < 2C. 
Applying an excitation above the band, with uj > 2C, it 
is possible to find the corresponding exact solutions for 
localized modes pinned to the VT complex: 



4>n 



VNja+iPb) ( x /w 2 +4C 2 +w S 



-M 



2uj - ^lu 2 + 4C 2 



2C 



(19) 



where ipA and ipB are given by the above solutions for 
the symmetric and asymmetric modes, i.e., severally, 
Eqs. p} and ([5), ([6]), with V and E replaced by 

V = V 2 / (2td - Vw 2 + \b 2 ) and E = E + V. This 
means that the solution for the nonpropagating symmet- 
ric modes remains explicit, while Eq. ^ for the asym- 
metric mode takes the form of a quartic equation for was • 
In conclusion, we have introduced the species of VT- 
symmetric systems with the balanced nonlinear gain and 
loss. For the dimer system, we have produced a complete 
set of analytical solutions, which feature the spontaneous 
symmetry breaking and multistability, not achievable 
in previously studied "PT-symmetric systems. We have 



demonstrated, also in the analytical form, that the sym- 
metric and asymmetric excitations in the dimer, if it is 
coupled to the linear chain, give rise to a variety of novel 
nonlinear Fano resonances, including the bistability be- 
tween them, and simultaneously to perfect-transmission 
regimes resembling EIT, as well as to the resonantly- 
amplified transmission. The coexistence of these scat- 
tering channels suggests an application to the design of 
data-processing schemes. Nonpropagating modes in the 
chain, pinned to the VT scatterer, were found too. 
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